Mon. Not. R. Astron. Soc. 000, ITHT21 (20091 Printed 11 February 2010 (MN M^jX style file v2.2) 



Irregular satellites of Jupiter: Capture configurations of 
binary-asteroids 

H. S. Gaspar*, O. C. Winter* and E. Vieira Neto* 

UNESP Univ Estadual Paulista, Grupo de Dinamica Orbital e Planetologia, CEP 12. 516-410, Guaratinguetd, SP - Brazil 



ABSTRACT 

The origins of irregular satellites of the giant planets are an important piece of the giant 
"puzzle" that is the theory of Solar System formation. It is well established that they 
are not in situ formation objects, around the planet, as are believed to be the regular 
ones. Then, the most plausible hypothesis to explain their origins is that they formed 
elsewhere and were captured by the planet. However, captures under restricted three- 
body problem dynamics have temporary feature, which makes necessary the action 
of an auxiliary capture mechanism. Nevertheless, there not exist one well established 
capture mechanism. In this work, we tried to understand which aspects of a binary- 
asteroid capture mechanism could favour the permanent capture of one member of 
a binary asteroid. We performed more than eight thousand numerical simulations of 
capture trajectories considering the four-body dynamical system Sun, Jupiter, Binary- 
asteroid. We restricted the problem to the circular planar prograde case, and time of 
integration to 10 4 years. With respect to the binary features, we noted that 1) tighter 
binaries are much more susceptible to produce permanent captures than the large 
separation-ones. We also found that 2) the permanent capture probability of the minor 
member of the binary is much more expressive than the major body permanent capture 
probability. On the other hand, among the aspects of capture-disruption process, 4) a 
pseudo eastern-quadrature was noted to be a very likely capture angular configuration 
at the instant of binary disruptions. In addition, we also found that the 5) capture 
probability is higher for binary asteroids which disrupt in an inferior-conjunction with 
Jupiter. These results show that the Sun plays a very important role on the capture 
dynamic of binary asteroids. 

Key words: planets and satellites: formation - minor planets, asteroids - Solar 
system: formation 



1 INTRODUCTION 

The existence of more than 350 natural satellites is known, 
from which approximately 50 % are planetary ones. An in- 
teresting point about this number, is that before 1997 just 
a tenth of such objects was known, i.e., the new "CCD ob- 



formation through matter accretion from the circumplane 



magnitude within iust a half decade (IGladman et al. 


200d. 120011: ISheppard & Jewrd 120031; 


Holman et al. 


Kavelaars et al.l 12004: Sheppard et al. 


20051. 2006ft. 



The planetary satellites can be distinguished into two char- 
acteristic g roups: regulars and irregulars (|Kuiperl 1 19561 : 
IPeald [ 1999). The first group, is characterized by small 
values of semi-major axis, eccentricities and inclinations. 
These characteristics are a strong signature of in situ- 
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tary d isc (ILunine fc Stevenson 
200l| : ICanup fc Ward! |2002l. 



1 19821 : IVieira Neto fc Winter 



2001 



6: iMosqueira fc Estrada 



20031 : ISheppard fc Jewittl \20o3 l In contrast, the irregu 



lar s atellites have large values of semi-major axis (I Burns! 
1 1986ft . often high eccentricities and inclinations. A large 
part of irregular satellites have r etrograde orbital incli- 
natio ns higher than 90 degrees (|jewitt fc Haghighipourl 
2JD07i). Another important characteristic of the irregular 
ones are the family groups, i.e., sa tellite groups charac- 
terized by similar or bital elements |GIadman et al.l 1200 ll : 
iKavelaars et al.ll2004T ). These peculiar characteristics are in- 
compatible wi th the in sit^- formation model through mat- 
ter accretion (|Kuiper| Il956l h and since they are the ma- 
jority group of planetary satellites in the solar system, 
there exists a large scientific interest about their origin. 
Then, the most plausible hypothesis to explain their ori- 
gins is that they formed elsewhere and were captured 
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by the planet dKuiper 19561; Heppenheimer fc Pored 1 19771 ; 
IPollack et ail 1 19791 ; IColombo fc Franklinl Il97lf l. However, 
many studies have shown that gravitation al captures un 
der three-body-dynami c s are t emporary dEverharj 
Heppenheimer fc Porco 1977 1; Carusi fc Valsecchi 



Benner fe Mckinnonl I 



Winter fc Vieira Netol 



991" 



1973 



1979 



200 ll ) 



IVieira Neto fc Winter! 12001 
This fact has induced re- 
searchers to propose some auxiliary capture mechanism. 
Among others, we point out four mostly well known: 



(i) Gas drag capture (|Pollack et all Il979l ; ICuk fc Burns! 

120041 '): A temporarily captured asteroid becomes perma- 
nently captured through kinetic energy decrease due to gas 

drag inside the circumplanet ary disk of gas and dust ; 

(ii) Pull-Down capture jHeppenheimer fc Pored 1 19771 ; 
IVieira Neto et aL | |2004| ; lOliveira et al.ll2007l ): A temporary 



captured asteroid becomes permanently captured due to an 
increase of the Hill's radius of the planet. This increase in 
Hill's radius occurs due to either planet's mass growth or 
planet's migration away from the Sun; 

(iii) Close-approach interaction captures 



(IColombo & Franklinll97ll;lTsuill999l, 20001; 


Astakhov et al.l 


2003; Ncsvornv et all 20031; Funato et all 


2004): A tem- 



porarily captured asteroid becomes permanently captured 
through energy and angular momentum exchanges with an 
existing satellite; 

(iv ) Capture of binary-ast eroids (|Agnor fc Hamilton! 
120061 : IVokrouhlickT et al.ll2008l ): One member of a binary- 
asteroid becomes permanently captured when the binary 
approaches the planet and disrupts. 

The capture mechanism of binary-asteroids is very in- 
teresting since the present observations have shown an in- 
creasing number of such systems in the main populations of 
such objec ts as the K u iper Belt, Main Bel t and N ear Earth 
Asteroids (iNolllbOOrj ). lAgnor fc Hamilton! l|2006h presented 
numerical simulations of close encounters between Neptune 
and a binary-asteroid, where they considered one asteroid 
comparable to Triton and a secondary, with equal mass or 
one order of magnitude lower. Their results show that is pos- 
sible to disrupt the binary when the close approach happens 
inside a spherical region whose radius, called tidal radius, is 
given by: 



Ttd = O.B 



( 3M P 

\mi + 7712 



1/3 



(1) 



where ap, Mp, mi and m,2 are the binary semi-major axis, 
planet mass, primary and secondary asteroid masses, respec- 
tively. A possible outcome after disruption is the capture of 
one member of the primordial binary-asteroid. 

The i ncreasing num ber of binary-asteroid dis coveries 
(|Nolll2006h along with the lAgnor fc Hamilton! (|2006h results, 
have motivated us to study the binary-asteroid capture pro- 
cess in the context of four-body dynamics, where we consid- 
ered Sun, Jupiter and a pair of asteroids. The purpose of this 
work is to identify the most important orbital characteristics 
inherent to the binary-asteroid capture/disruption process 
that produce the permanent capt ure of at least one mem- 
ber. Compared to existin g works (jAgnor fc Hamiltonll2006l : 
IVokrouhlickv et alj |2008h. our study considers the inclusion 
of solar perturbation. We found that the Sun's presence has 
a crucial influence on the binary capture/disruption process, 
at least in the planar case. 



This paper is built with the following structure: Section 
[2]describes the model we use in our study as well the adopted 
numerical approach. Section [3] presents the results with an 
analysis of them. Finally, the last section summarizes our 
conclusions. 



2 CAPTURE MODEL 

Given the temporary feature of captures in the three-body 
problem, we have chosen to perform a study under four-body 
dynamics using the Sun and Jupiter as primary bodies and 
a binary-asteroid. We basically propose a capture model in 
which a binary-asteroid first becomes temporarily captured 
by Jupiter, and then disrupts and has one of its member per- 
manently captured by Jupiter while the other one escapes. 
The present paper addresses the early results of a more gen- 
eral work which is under development. It should be noted 
that the main goal of the present work is not to reproduce 
the actual configuration of Jupiter's irregular satellites, but 
rather, to comprehend how specific configurations can lead 
an asteroid member of a primordial binary to a permanent 
capture. As a first stage, we have considered only the pla- 
nar prograde case. Given this capture model, our task is to 
search for the initial conditions which yields asteroid tempo- 
rary captures by Jupiter. The temporary captures, as well 
as the close encount ers, are intrins i c features of solar system 
formation theories (Pollack ct al. 1996; H ahn fc Malhotra 
l2005l ; lTsiganis et al.ll2005l ; lGomes et alj|2005h . Furthermore, 
temporary captures seem to be a more efficient mechanism 
to accomplish asteroid captures because there is a longer in- 
teraction time between the binary-asteroid and the planet 
rather than a single passage with a very short time of inter- 
action. 



2.1 Procedure 

In this work, as a first study, we considered only the copla- 
nar four-body dynamics. Furthermore, we set Jupiter's ec- 
centricity to zero for all the simulations. In order to perform 
the numeric stu dies we used an integrator based on Gauss- 
Radau spacing (|Everhartll 19851 ). In order to verify the inte- 
grator's accuracy we checked whether the system's total en- 
ergy holds throughout the integration. We found the energy 
variation was lower than 10 -11 . In addition, we have checked 
the value of the Jacobi constant of the binary-asteroid cen- 
ter of mass before binary disruptions for all the trajectories. 
We found that the Jacobi constant variation was lower than 

10 -9 

The adopted procedure basically followed three phases: 
i) Firstly, we performed a capture time analysis of the sys- 
tem Sun- Jupiter-particle through which we obtain the, from 
now on designate, suitable initial conditions, i.e., initial con- 
ditions which result in a particle's temporary capture by 
Jupiter, ii) Given the suitable initial conditions, we replace 
the individual particle by a pair of bodies in order to set 
up a binary-asteroid, which will be called initial conditions. 
iii) Finally, we study the binary-asteroid capture through 
simulations of the system Sun, Jupiter, binary-asteroid, by 
considering a set of different initial conditions derived from 
each one of the suitable initial conditions. 
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Figure 1. Example of escape trajectory in backward integration. 
The particle is in clockwise trajectory under Jupiter fixed frame 
of referenceThe initial point indicated by blue 'x' corresponds to 
the primary initial conditions (around Jupiter). The final point 
indicated by red '+' corresponds to the suitable initial condition, 
which will be replaced by the binary-asteroid. 



2.2 Suitable Initial Conditions 

In order to obtain the suitable initial conditions, we per- 
fo rmed a primary captu r e tim e study following the steps 
of IVieira Neto fc Winter! (120011 ). It consists of the integra- 
tion of a particle trajectory using a negative time step un- 
der the three-body dynamics considering Sun and Jupiter as 
primaries. By setting the particle to start orbiting around 
Jupiter we have three possible outcomes depending on par- 
ticle's initial condition: i) The particle collides with Jupiter, 
ii) The particle remains orbiting the vicinity of Jupiter up 
to the final time of integration (10 4 years), in such cases, the 
particle's initial conditions are stable ones, iii) The particle 
escapes from Jupiter's vicinity and begins to orbit the Sun. 

In order to make easy the comprehension, we finally 
define these particle's orbital elements around Jupiter as 
primary initial conditions. Summarizing: 

Primary initial conditions: Three-body problem initial 
conditions which are backward integrated in time in order 
to obtain the suitable initial conditions; 

Suitable initial conditions: Three-body problem initial 
conditions which lead the particle to temporary captures 
by Jupiter. By replacing the particle for a binary-asteroid 
one derives the initial conditions of the binary-asteroids; 

Initial conditions: The real initial conditions used to study 
the binary-asteroid capture dynamics, in which we consider 
Sun, Jupiter and a pair of asteroids; 

We are particularly interested on the data relative to escape 
trajectories, which are capture ones when it is considered 
time forward. The particle's and Jupiter's orbital elements 
with respect to the Sun at the instant when the particle 
escape^J are taken as the suitable initial conditions, as illus- 
trates the example in Figure [T] 

The main results of our three-body dynamics simula- 
tions, Sun- Jupiter-particle, is shown in the capture time map 



1 We check the particle's two-body energy every time interval of 
one year during the numerical simulations. We consider that the 
escape occurs when the particle's two-body energy with respect 
to Jupiter becomes positive. 
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Figure 2. (a) Capture time mapping. The color scale correspond 
to the capture time (in years) for such pair of initial semi-major 
axis and eccentricity. The green '+' indicates the initial conditions 
which produce collisions. Initials longitude of pcricenter and true 
anomaly of the particle were set to zero (ro = / = 0). (b) Primary 
initial conditions used to obtain the suitable initial conditions. 
The green and black 'x' indicate the cases whose capture time 
are shorter and longer than a thousand of years, respectively. 
The pink dotted line is given by Equation [2] 

of Figure [2ja). This is an axe diagram whose color scale de- 
notes particle's escape time (in years) in the backward's in- 
tegration. In this plot, the particle's semi-major axis is given 
in terms of Hill's radius of Jupiter. This map has been gen- 
erated by setting the particle's initial longitude of pericenter 
(w) and initial true anomaly (/) with respect to Jupiter as 
zero. 

In plot (a) of Figure [2] the yellow region corresponds to 
the primary initial conditions which do not result in escapes 
from Jupiter in 10 4 years. Considerin g the satellite eccentric- 
ity up to 0.5. lDomingos et al.l (|2006t ) obtained an expression 
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for this stable region as a critical semi-major axis inside 
which the satellites would remain stable. This expression is 
given by 



0.4895(1.0000 - 1.0305e P - 0.2738e sa t) 



(2) 



where ep and e sa t are the planet's and satellite's eccentrici- 
ties, respectively. Particularly, in this work, the second term 
between parenthesis on the right hand side of Equation [2] 
vanishes si nce we have taken ep = 0. In order to obtain this 
expression IDomingos et all (120061) fo llowed the same proce- 
dure used by Vieira Neto fe Winter! |200ll ) and also set the 
particle's initial longitude of pericenter (zn) and initial true 
anomaly (/) as zero. In section [3.31 we will show that the 
stable region can be more extensive for initial values of vj 
and / different from zero. 

Finally, the non-yellow region on the map (a) of Fig- 
ure [2] indicates the primary initial conditions that resulted 
in escape. These are the cases from which we can take the 
suitable initial conditions. However, it is not feasible to use 
all the data obtained from the capture time analysis. By an- 
alyzing the capture times of escape cases, we found only 45 
cases in which the capture time exceed one thousand years. 
It is feasible to take all these long time capture cases to com- 
pose the set of suitable initial conditions. Among the cases 
in which the capture times are shorter than a thousand of 
years, we took an uniformly-spaced grid of points in a x e 
space in order to complete the set of suitable initial condi- 
tions with a representative set of the whole data. Plot (b) 
in Figure [2] summarizes the set of primary initial conditions 
we used to obtain the suitable initial conditions. The axe 
diagram of Figure [3] shows the suitable initial conditions 
which was obtained from the primary ones. The plot shows 
the semi-major axis and eccentricities of the main asteroid 
in the heliocentric frame. Furthermore, all the trajectories 
are direct with respect to the Sun. It is also plotted, in Fig- 
ure |31 the Tisserand rel ations for T = 2 .996 and T = 3.036. 
The Tisserand relat i on |Tisserandlll896T l. as one can find in 
iMurrav fc DermottJ (|l999T l. is given by: 



e 2 ) cos(7) ss constant 



(3) 



where a, e and / are the object heliocentric semi-major axis, 
eccentricity and Inclination. 



2.3 Model for the Binary Asteroid 

The term binary-asteroid 

can refer to either a system in which a pair of asteroids of 
similar masses orbit their common b arycenter or an asteroid 
that has a small satellite (|Nonll2006r i . In order to compose a 
binary-asteroid by using the suitable initial conditions, we 
added a second body orbiting the first one. From now on, the 
main and the secondary asteroids will be called PI and P2. 
For simplicity, we model a binary-asteroid by setting P2's 
orbital elements with respect to PI. In other words, asteroid 
P2 initially orbits asteroid PI. P2's orbital elements with 
respect to PI will be referred as binary's elements. Given 
our particular interest in Jupiter's irregular satellites, we 
set the P2's mass as mi — 10 19 kg based on Himalia's mass , 
the largest irregular satellite of Jupiter ([Em clvaiiov 2005h . 
Using the same mass ratio as lAgnor fc Hamilton! ( 20061 ) . we 
set Pi's mass mi = 10rri2 = 10 20 kg. 
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Figure 3. Suitable initial conditions axe heliocentric diagram. 
The green "x" and black "+" indicate suitable initial conditions 
which were obtained from primary initial conditions whose cap- 
ture time are shorter and longer than a thousand years, respec- 
tively. Red and blue curves are Tisserand relations T = 2.996 and 
T = 3.036, respectively, which encompass all the suitable initial 
conditions. 



As already stated, the main goal of this paper is to iden- 
tify the most appropriate configurations that would gener- 
ate permanent capture of one asteroid from a binary system. 
Thus, since we have a huge range of possibilities, is this pa- 
per we made some restrictions in order to be able to explore 
a significant part of the initial conditions space. Among the 
restrictions, we considered that P2 is always initially in pro- 
grade circular orbit around PI, (e = I = 0). 

From each of the 81 suitable initial conditions we de- 
rived 108 new initial conditions. We vary the initial binary's 
true anomaly /p, from up to 330° in steps of 30°, and 
the initial binary's semi- major axis ap, from 0.1 9~h up to 
0.5 rh in steps of 0.05 tv Here r^, with 'h' written in low- 
ercase, is Pi's Hill's radius with respect to the Sun calcu- 
lated for each one of the suitable initial conditions. Let's 
make clear that we did not make use of the tidal radius 
given by Equation [ T]The Hill's radius, as one can find in 
IMurrav fc Dermottl (|l999h . is defined by: 



(!) 



1/3 



(4) 



where fi and a are the mass ratio and the semi-major axis, 
respectively. The chosen upper semi-major axis limit of 0.5 
r~h is a well-established limit of stability for prograde sy stems 
(IHamilton fc Burnslll99ll . Il992l : IDomingos et al.ll2006r i . The 
lower semi-major axis limit was arbitrarily chosen, though 
tighter binaries do exist. By calculating the Hill's radius of 
PI Th with respect to the Sun (Equation 2}, one finds that 
these values vary from 1 x W~ S AU to 2 x 10~ 3 A[/., and P2's 
initial orbital velocity with respect to PI vary from 7m /s to 
22m/s. 



2.4 Binary asteroid's capture simulations 

Since we derived 108 initial conditions from each one of 81 
suitable initial conditions, we performed a total of 8 748 
binary-asteroid capture trajectory simulations. We set the 
trajectory integration time to 10 4 years and the output time 
step to 10 hours. In order to avoid a large amount of data 
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storage we designed an algorithm which identified the inte- 
gration main stages. We labeled each instant of these main 
stages, as follow: 

T\\ instant when the binary-asteroid is first captured by 
Jupiter; 

T2: instant when the binary-asteroid disrupts; 
T3: instant when only one member of the disrupted 
binary-asteroid escapes from Jupiter; 

Once the algorithm identifies each one of the three in- 
stants, it stores the instantaneous system configuration data, 
as well as the integration instant, in three distinct files. By 
taking the difference between Ti and either T3 or the fi- 
nal time of integration, we can compute the capture time 
for each case. This identification algorithm basically con- 
sists on a two-body energy check-up, every integration step, 
described as follow: 

The binary-asteroid approaches Jupiter in a quasi- 
parabolic trajectory, given that it initially orbits the Sun. 
For all the cases considered in this study at least one aster- 
oid two-body energy with respect to Jupiter was positive. At 
the instant we found that both Pi's and P2's two-body en- 
ergies with respect to Jupiter became negative the instant 
Ti is identified. Although, those energies do not become 
negative at the same instant. 

Similarly, P2's two-body energy with respect to PI is 
initially negative given that P2 initially orbits PI. Thus, if 
P2's two-body energy with respect to PI becomes positive, 
instant T2 is identified. Otherwise we do not identify neither 
instant T2 nor T3. Rather it could happen a double capture, 
a mutual collision or a double escape. However, in all the 
cases considered in this study the binary disrupted or collide 
with each other. 

Therefore, after the binary-asteroid's capture by 
Jupiter, both Pi's and P2's individual two-body energies 
with respect to Jupiter are negative. Furthermore, after the 
binary-asteroid's disruption, interactions between PI and P2 
become negligible, allowing either PI or P2 to individually 
escape from Jupiter. Finally, at the instant in which either 
Pi's or P2's two-body energy with respect to Jupiter became 
positive the instant T3 is identified. 

Three possible outcomes succeed T3: 

(i) The remaining asteroid collides with Jupiter, which 
characterizes a collision; 

(ii) The remaining asteroid escapes from Jupiter, which 
characterizes a double escape; 

(iii) The remaining asteroid stays bound throughout the 
rest of the integration time, which we characterize as a per- 
manent capture; 



3 RESULTS 

We show in this section some plots built with the data stored 
at instants Ti, T2 and T3 in three respective subsections. 
We will discuss some statistical results in the fourth subsec- 
tion, and show some examples of capture trajectories in the 
last subsection. 
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Figure 4. Binary's initial semi-major axis versus variation of 
the binary's semi-major axis at instant T\ for the cases which 
resulted in permanent captures of either PI or P2, in red and 
blue triangles, respectively. Red and blue filled lines are Pi's and 
P2's permanent capture histograms, respectively. 



3.1 Analysis at the instant of binary capture (Ti) 

The plot in Figure[4]compares the binary-asteroid separation 
at instant Ti with its initial separation, through binary's 
semi-major axis aB at instant T\ and initial binary's semi- 
major axis Ob, respectively. Both axis in this plot are mea- 
sured in units of initial Hill's radius of PI rh. This plot al- 
lows us to comprehend the binary-asteroid's evolution from 
the beginning of the integration to instant T\. 

The histograms reveal that i) tighter binary-asteroids 
are more susceptible to permanent captures than binary- 
asteroids with large separations and that ii) the capture 
probability of the binary's smaller member is much higher 
than the major companion's capture probability. The blue 
histogram shows a roughly negative exponential behavior 
of permanent capture probability, with respect to the ini- 
tial binary separation, which becomes very low for initial 
as ~ 0.35 TV Furthermore, beyond initial ob = 0.35 7~h one 
can observe that the permanent capture occurs more of- 
ten with binaries whose semi-major axis decreased. Con- 
sequently, these results point out to a limit as ~ 0.4 r^ be- 
yond which permanent capture plausibility is negligible. Ev- 
idently, Aob dispersion increases as the initial separations 
increases since it causes the binary bound to be weaker, con- 
sequently more susceptible to secular variations due to Sun 
and Jupiter perturbations. 

As weaker bounded binaries disrupt more easily, one 
could expect that they would more easily generate perma- 
nent captures. However, our results show that tighter bina- 
ries have higher probability to generate permanent captures. 
This apparent paradox can be understood in terms of the 
energy exchange needed to turn a temporary captur e into a 
permanent one. Based on results bv lTsuil (1 19991 . |2000| ). which 
show that it is possible for an asteroid to be kept captured 
by a planet due to exchange reactions with a local satel- 
lite, it is possible to explain the binary capture mechanism 
through energy exchanges in four steps: 

(i) Lets consider a binary-asteroid, initially orbiting the 
Sun, which will be captured by Jupiter. Since the binary- 
asteroid is primordially orbiting the Sun, each one of its 
members' individual energy is higher than the escape energy 
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so, i.e., the minimum necessary energy to allow each asteroid 
individually escape from Jupiter. 

(ii) However, once the two asteroids orbits closely their 
common barycenter, angular momentum and energy ex- 
changes occurs constantly. Furthermore, after the binary- 
asteroid be temporarily captured (Ti), Jupiter starts to dis- 
turb the binary binding. Therefore, the exchange reactions 
become more intense; 

(hi) As a consequence of Jupiter's perturbation the 
binary-asteroid disrupts (T2). However, before the binary 
disruption the energy exchanges between the asteroids pro- 
vides some energy states in which one member's energy is 
lower than the escape energy eo, which will not allow it to 
escape from Jupiter; 

(iv) Finally, after the binary disruption the interactions 
between the asteroids become negligible, so that, the aster- 
oid whose energy decreased to values lower than eo remains 
captured by Jupiter while the other whose energy increased 
will escape from Jupiter after some time (T3). 

Note that these ideas agree well with our results: 

(i) The rupture of tighter binaries imply on larger energy 
exchange, implying on higher probability of permanent cap- 
ture; 

(ii) The smaller body of the binary is the one that suffers 
larger energy exchange and consequently is the one with 
higher probability of permanent capture; 

(iii) Finally, the the existence of a separation limit agrees 
well with the conclusions since the binary separation is pro- 
portional to the binding energy, which, in its turn, corre- 
sponds to the maximum energy that the asteroids can ex- 
change. In other words, larger separation-binaries can not 
provide enough energy exchange between its members in or- 
der to allow one of them to became permanently captured. 

Plot (b) of Figure[S]illustrates this energy exchange pro- 
cess for the trajectory shown in plot (a). Plot (c) shows the 
time evolution of the Jacobi constant value for each individ- 
ual asteroid. It illustrates that shortly after instant T2 the 
interaction between PI and P2 becomes negligible. Note also 
that, the Jacobi constant value of P2 higher than Jacobi con- 
stant value of Lagrangian point L\. Therefore, shows that 
P2 will never escape Jupiter's vicinity 

Nevertheless, we could expect this increasing fraction of 
captured low-semimajor axis binaries to have a maximum 
where it must roll over since very tighter binaries should 
reach a limit in which the binary-asteroids can be consid- 
ered as a particle and would never disrupt. Furthermore, we 
should also expect an increasing fraction of mutual collisions 
inasmuch binary semi-major axis decreases. 

3.2 Analysis at the instant of binary disruption 

(Ta) 

As defined before, the instant T2 is characterized by the 
binary-asteroid disruption. So, the graphics in this subsec- 
tion refer to elements of each asteroid individually with re- 
spect to Jupiter. 

Plot (a) in the Figure|6]is an axe diagram, at the instant 
T2 , of the asteroids that remained permanently captured by 
Jupiter. Plot (b) shows the axe diagram of the last asteroid 
to escape from Jupiter for double escape cases. By compar- 



ing these two diagrams, one finds a region on axe space in 
which an asteroid remains permanently captured if its semi- 
major axis and eccentricity are enclosed within it at instant 
T2. It must be clear that Figure [6] shows T2 instantaneous 
diagram of osculating elements which must vary in time due 
to solar perturbation. Furthermore, we should expect some 
weak interaction between the pair so far as they recede suf- 
ficiently away from each other. By checking the variation of 
Jacobi constant value, of each individual asteroid, we could 
estimate how long it takes to this mutual interaction become 
negligible. By considering a variation of the order of 1CF 8 
in Jacobi constant value, we found that for the worst case it 
took about half orbital period about Jupiter (~ 170 days) 
to satisfy the condition. 

By fitting an expression that bounds this region in the 
axe space we found a limit on semi-major axis given in 
terms of eccentricity, as follow: 

a*(e) = 0.4500(1.0000 - 0.2046e) (5) 

The condition a < a*(e) at instant T2 can be thought 
as sufficient but not necessary capture condition. In other 
words, if the captured object obeys a < a*(e), then our sim- 
ulations show that the temporary capture always becomes 
permanent. However, we also found large number of perma- 
nent captures that had a slightly greater semi-major axis in 
a region where there is a mix of captures and double escapes 

A second graphic at instant T2 is as an angular illus- 
trative histogram. One can better understand the angular 
configuration analyzing the phase angles 9i and 62 shown 
in Figure [7] The angle 8\ gives the PI phase angle with re- 
spect to the Sun- Jupiter direction, while, 62 gives P2 phase 
angle with respect to the Jupiter-Pl directions. Note that 
in this system PI orbits counterclockwise Jupiter and P2 
orbits counterclockwise PI. 

Analysing the angular histograms of Figure [5] we see 
that: i) the disruption preferentially occurs when Jupiter, PI 
and P2 are approximately aligned, i.e., 62 « or 62 ~ 180°. 
The smaller probability for 81 ~ 180° in capture histogram 
(a) , plus the higher probability observed for 9\ « 180° in es- 
cape histogram (b) indicate that ii) disruptions which occurs 
when binary-asteroid is located aligned between Jupiter and 
Sun most likely result in double escapes. On the other hand, 
histogram (a) shows that iii) permanent capture of P2 most 
likely results from binary-asteroids which disrupt at a angu- 
lar position approximately 90° after it cross the Sun-Jupiter 
line (9i ~ 270°). Finally, from the capture histogram (a) 
one also notes that permanent captures of P2 succeed from 
disruptions which occurred when P2 were at inferior con- 
junction with PI, as seen from Jupiter (82 ~ 180°). 

The preference for captures with the smaller asteroid 
unbinding when closer to Jupiter (82 180°) is simply un- 
derstood as being due to the velocity vector about its center- 
of-mass motion being opposite to the planetocentric velocity 
in that geometry. Consequently, asteroid P2 gets its individ- 
ual speed, with respect to Jupiter, reduced to a minimum. 

The preference for captures when 8\ ~ 270° is under- 
stood as being due to the solar perturbation, i.e., the Sun's 
gravity tends to increase the velocity of binary center of 
mass about Jupiter when 61 » 90° while it tends to de- 
crease this velocity when 61 ~ 270°. These characteristic 
angular positions with respect to the Sun- Jupiter line, rein- 
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Figure 5. A first example of capture process. In plot (a), red '+' and blue 'x' show the trajectories of PI and P2 in Jupiter's (black 
circle) planctocentric non-rotating frame of reference, respectively. The output frequency is 10 days. Dashed grey line and the black 
arrow, in the zoom window, indicate Sun direction at the disruption instant (T2). Tracking binary-asteroid trajectory, one sees that it 
becomes temporarily captured (light blue triangle) by Jupiter, disrupts (green square inside zoom box) and finally has its minor member 
permanently captured by Jupiter while its major member escapes from Jupiter's vicinity (pink down triangle). Plots (b) and (c) show 
the time evolution of energy and Cj for the trajectory of plot (a), respectively. In plot (b), red '+' and blue 'x' are the two-body energies 
of PI and P2 with respect to Jupiter, respectively, and the pink filled line is the binary's two-body energy, i.e., P2's two-body energy 
with respect to PI. In plot (c), red '+' is the difference of Cj values calculated for asteroid PI and L\ Lagrangian point, and similarly , 
blue 'x' difference for P2 and L\. 
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Figure 7. Sketch of angular configuration between the bodies. 
One can infers a Sun-Jupiter-Pl alignment if 8i=0 or #i=180°and 
a Jupiter-Pl-P2 alignment when 6*2=0 or 02=180°. 
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Figure 6. Diagrams of semi-major axis versus eccentricity of the 
asteroids at the instant T2. Plot (a) shows the cases which re- 
sulted in permanent captures of either PI or P2, in red and blue, 
respectively. Plot (b) shows the cases which resulted in double 
escapes. The green '+' corresponds to the orbital elements of the 
asteroids that escape after instant T3. Pink dotted straight line 
is the fitted limit semi-major axis given by Equation [5] 



force the importance of the solar presence on the dynamics 
of binary-asteroid captures: 



3.3 



Analysis at the escape instant of one asteroid 
(Ts) 



The data stored at instant T3 shows the final configuration 
of the captured asteroid, given that the captured asteroid 
will not suffer any interaction with its primordial partner. 
Figure [9] is an a x e diagram of captured asteroid at instant 
T3. The pink filled contour is an extended border of a more 
general region of stability. In fact, we firstly worried about 
the points located bey ond critical semi-major axis found by 
iDomingos et al.l (|2006l ) since they represent the permanently 
captured asteroids. Nevertheless, by performing a more gen- 
eral capture time analysis we found the extended stability 
border. 

We performed this more general study following the 
same procedure described at section [2T2l but considering dis- 
tinct initial values for longitude of pericenter voq and true 
anomaly fo. From the results, shown in the Figure [lOl we 
found that the pairs of initial values (zuq — 0, fo = 180°) 
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Figure 8. Angular histograms of instant T%, The red circles 
indicate the angular position of PI with respect to Sun-Jupiter 
direction given by angle 6\. The angular sections that surround 
each red circle indicate the angular position of P2 with respect to 
Jupitcr-Pl direction, given by 82 (see Figure [71. The color of the 
angular sections in histogram (a) correspond to the percentages of 
the total number of P2's permanent captures (972; see TableflJ, 
and the color of angular sections in histogram (b) correspond to 
the cases which resulted in double escapes (7446; TableflJ, given 
by the respective color scales. White color indicates absolute null 
number of events. The bodies orbit counterclockwise. 
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Table 1. Percentages of captures and collisions. 



Description 


Short capture times 


Long capture times 


Simulations 


3 888 


4 860 


Captures 


7 


972 


Captures of PI 


1 


5 


Captures of P2 


6 


966 


Double captures' 1 





1 


Collisions 


180 


113 



a PI and P2 captured 

and (wo = 180°, /o = 180°) yield capture time maps in 
which the stability regions are much larger than for initial 
values (vjo = 0, fo = 0), shown in Figure [2] By combining 
the borders of both maps, in such manner we obtain the 
larger region of stability, we built the referred more general 
stability edge shown, as a pink filled line, in Figure [5] 

The points in the plot of Figure [5] shows that the final 
orbital elements of permanently captured asteroids cover a 
wide region in the axe space. Most of the captured objects 
are very far from the planet (ae > 0.4 th) and will probably 
be removed due to perturbations not included in our study. 
The good candidates to survive are those closer to the planet 
(ob < 0.35 m). Since in this study we considered only the 
prograde planar case, we are not able to compare our results 
with the currently known Jupiter's irregular satellites. How- 
ever, we included them in the plot of Figure just to have 
an idea of the orbital shape of captured objects. In order to 
make a fair comparison it will be needed to make a study in 
the 3-D space considering the inclinations. This is a study 
in progress 

3.4 Binary captures in numbers 

Among 8 748 simulated trajectories, 4 860 are from long cap- 
ture time primary initial conditions and 3 888 from short 
capture time ones. Table [T] presents the quantities of perma- 
nent captures and collisions. Second and third columns show 
the values with respect to the short and long capture time 
cases, respectively. As it is shown in Table[T] though perma- 
nent capture probability of the cases derived from short-time 
primary initial conditions are low (0.18%), the permanent 
capture probability of cases derived from long time condi- 
tions are much larger (20%). The collision probabilities for 
both long and short times derived from primary initial con- 
ditions have the same order and are not negligible. 

3.5 Sample of capture trajectories 

This section presents some examples of capture trajecto- 
ries from three distinct cases. Firstly, we already showed in 
Figure [^a) , a typical example where a binary-asteroid ap- 
proached Jupiter, became captured, disrupted after a while 
and had its minor member permanently captured by Jupiter 
while its major member escaped Jupiter's vicinity. Fig- 
ure [TTJ a), shows rare example where the major asteroid re- 
mains captured by Jupiter after disruption of the primordial 
binary-asteroid. Furthermore, another peculiarity in this ex- 
ample is T2—T3, i.e., the escape of P2 happens at the in- 



stant of binary disruption. Finally, Figure fTTT b) presents our 
single case among > 8000 performed where both asteroids 
remained captured by Jupiter up to the end of integration, 
even after the disruption of the primordial binary-asteroid. 
In both examples of the Figure [11] we simulated the tra- 
jectories for 10 4 years, and for both cases the instant T2 
happens before 10 2 years of integration. Nevertheless, by 
checking the Jacobi constant value for each one of the bod- 
ies of the double capture example (Figure ITTb ). we found 
that their values are smaller than the L\ Jacobi constant 
value. Therefore, these bodies may eventually escape from 
the planet. 



4 CONCLUSIONS 

In this work we have studied the capture dynamics of binary- 
asteroids by looking for favorable conditions of capture. 
Our results present new perspectives about the problem 
of binary-asteroid captures emphasizing the importance of 
Sun's role in the dynamics. The Sun is not necessary to pro- 
duce a binary rupture since Jupiter alone can do it. However, 
the Sun plays a key role in the disruption process in order to 
make the capture to become permanent. The results allow 
us to comprehend about both binary-asteroid's features as 
well as intrinsic features of capture process' main stages. 

The observed characteristics at the first main stage, Ti, 
have revealed that: i) Tighter binary-asteroids are more sus- 
ceptible to permanent captures than binaries with larger 
separation. In fact, the permanent capture probability be- 
haves inversely proportional to the binary's semi-major axis. 
This results indicate that binary's energy exchange allows 
the asteroid to become permanently captured. That is, since 
tighter binaries interactions are more intense, its members 
can exchange higher amount of energy. In such a way one 
asteroid can have its energy sufficiently decreased in order to 
not be able to escape from Jupiter. Through this conclusion, 
we could argue that binary-asteroids with high eccentricities 
would disrupt more easily, but would not exchange the nec- 
essary amount of energy to result in a permanent capture. As 
mentioned, there must exist a lower binary-separation limit 
below which the binary never disrupts and consequently 
there is no capture. 

The observed characteristics at the second main stage, 
T2, tell about process' features. It was shown that the an- 
gular position of bodies at disruption instant (T2) are re- 
lated with the permanent capture probability. Summarizing: 
ii) Disruption preferentially occurs when both asteroids are 
approximately aligned with Jupiter; Nevertheless, iii) dis- 
ruptions which occurs when P2 is located between Jupiter 
and PI result more often in permanent capture; Finally, we 
found that iv) the permanent capture probability is higher 
when the binary-asteroid disrupts in an eastern quadrature 
i.e., an angular position approximately 90° after the binary 
cross the Sun- Jupiter direction. 

The permanent capture probability for an specific set 
of initial conditions, such derived from long time primary 
initial conditions, was shown to reach 20%. Therefore, the 
good candidates are those derived from long capture time 
primary initial conditions. 

As a note, we give here a reference of a similar work sub- 
mitted to Icarus journal, which also considers the solar per- 
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Figure 9. Same as plot (a) in Figure \E\ for instant T3, though. Red '+' and blue 'x' represent the cases which resulted in permanent 
captures of either PI or P2, respectively. Black circumferences represent the orbital elements of Jupiter's real prograde irregular satellites. 
The pink dotted line is the same as in plot (a) of Figure [6] The filled pink contour is a stability edge we found numerically as discussed 
in section [3.31 
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turba tion. This work is available on astro-ph ijPhilpott et al.l 
l2009h . 

Finally, as the main goal of this paper was to address 
the favorable conditions which makes the permanent cap- 
ture plausible, we have chosen a procedure to get initial 
conditions without taking into account where the incoming 



objects came from. It means that, in this work we have just 
tried the model plausibility without taking into account how 
it could reproduce the currently observed objects. So in or- 
der to get a more realistic probability, it must be considered 
several aspects as inclinations, mass ratios, binary eccentric- 
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Figure 11. Two examples of capture pathways. Red '+' and blue 'x' are the coordinates of PI and P2, respectively. Light blue triangle 
and green square denote Tiand T2, respectively, (a) The asteroid PI remains permanently captured by Jupiter, (b) Both asteroids 
remain captured by Jupiter even after disruption. 



ities as well as to study how realistic are the trajectories of 
incoming objects. 
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